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ABSTRACT 



The 'Bloated Stars Scenario' proposes that AGN broad line emission originates in the 
winds or envelopes of bloated stars (BS). Its main advantage over BLR cloud models 
is the gravitational confinement of the gas and its major difficulty the large estimated 
\l number of BSs and the resulting high collisional and evolutionary mass loss rates. Pre- 

vious work on this model did not include full calculations of the spectrum and used very 



^ simplified stellar distribution functions. Here we calculate the emission line spectrum 



by applying a detailed numerical photoionization code to the wind and by assuming 
a detailed nucleus model (Murphy, Cohn & Durisen 1991). We study a wide range of 
wind structures for a QSO model with L lon = 7xl0 45 erg/s, Mbh = 8xlO 7 M and 
^T) a stellar core density of 7xlO 7 M0/pc 3 . The BSs emit a spherically symmetric wind 

whose size and boundary density are determined by various processes: Comptonization 
by the central continuum source, calculated here self consistently, tidal disruption by 
the black hole and the limit set by the wind's finite mass. We find that the emission 
^~- ) spectrum is mainly determined by the conditions at the boundary of the line emitting 

fraction of the wind rather than by its internal structure. Comptonization results in 
f^S a very high ionization parameter at the boundary which produces an excess of unob- 

served broad high excitation forbidden lines. The finite mass constraint limits the wind 
ON size, increases the boundary density and thus improves the results. Slow, decelerating, 

mass-constrained flows with high gas densities (10 8 to 10 12 cm -3 ) are as successful as 
| cloud models in reproducing the overall observed line spectrum. The MgllA2798 and 

I NvA1240 lines are however under-produced in our models. We adjust the number of 

O BSs so as to obtain the observed EW(La) and find that only ~ 5 x 10 4 BSs with dense 

winds are required in the inner 1/3 pc. This small fraction approaches that of super- 
giants (SG) in the solar neighborhood and the calculated mass loss is consistent with 
the observational constraints. The required number of BSs, and their mass loss rate, 
are very sensitive functions of the wind density structure. SG-like BSs are ruled out 
by their high mass loss. BSs with dense winds can reproduce the BLR line spectrum 
and be supported by the stellar population without excessive mass loss and collisional 
destruction rates. The question whether such hitherto unobserved stars actually exist 
in the BLR remains open. 
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1 INTRODUCTION 

The observed properties of active galactic nuclei (AGN) lead 
to the conclusion that the line emission originates in numer- 
ous small, cold and dense gas concentrations, denoted by the 
general term clouds. It is inferred that the Broad Line Region 
(BLR) cloud system typically consists of ~ 10 r clouds with 
particle densities of 10 8 — 10 12 cm -3 and an ionized fraction 
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of dimension 3xl0 12 cm at a temperature of ~ 10 4 K. How- 
ever, the clouds' distribution and velocity field, as well as 
the functional dependence of a cloud's density and size on 
the distance from the black hole, cannot be deduced from 
present day observations and are therefore free parameters 
of the cloud model. Once these parameters are specified by 
additional assumptions, the BLR line emission can be calcu- 
lated by photoionization codes. These calculations demon- 
strate that such cloud models can reproduce to a fair degree 
the observed line properties (Netzer 1990). 
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The major objection to the clouds model is the con- 
finement problem. The clouds, with estimated masses of 
10 -9 Mq, are not self gravitating objects and should disinte- 
grate on the sound crossing time scale of ~ 0.5 yrs, whereas 
the crossing time of the BLR is ~ 30i^g 2 yrs (where L±e 
is the ionizing continuum luminosity in units of 10 46 erg/s). 
Unless clouds are continuously formed in the BLR, an ex- 
ternal confinement mechanism is necessary. A once promis- 
ing idea, the two-phase model (Krolik, McKee & Tarter 
1981) whereby pressure equilibrium is maintained between 
the cloud (cold phase at ~ 10 4 K) and the hot interstellar 
medium (HIM, hot phase at ~ 10 8 K) has met serious objec- 
tions. The gas temperature in such models depends on the 
shape of the UV-X-ray continuum. Realistic continua mod- 
els predict the HIM temperature to be ~ 10 r K which is not 
compatible with two-phase equilibrium (Fabian et al 1986; 
Mathews & Ferland 1987). Additional problems are cloud 
breakup by drag forces due to the supersonic motion of the 
clouds, the unobserved opacity expected from this cooler 
HIM and the instability of the cold phase due to internal 
line pressure. 

An alternative model that has not been investigated 
in detail yet is the 'stellar winds scenario' (Edwards 1980; 
Mathews 1983; Scoville & Norman 1988; Penston 1988; 
Kazanas 1989). This scenario proposes that the lines are 
emitted from the mass-loss winds or envelopes of super- 
giants, and is motivated by the following arguments: 1) The 
confinement problem is resolved by introducing a self grav- 
itating gas reservoir — the star — that can continuously re- 
plenish the cold gas that is evaporating to the HIM. 2) Stars 
are expected to exist near the central source and estimates 
suggest that the conditions in supergiant winds resemble 
those inferred for the line emitting gas. 3) The lack of ev- 
idence for radial velocity of the BLR gas (e.g. Maoz et al. 
1991) is consistent with the isotropic velocity field of the 
virialized stellar core around a massive black hole. 

Earlier work on stellar wind models was confined to 
simplified approximations. In particular, there are yet no 
detailed calculations of the line ratios and the line profiles 
were approximated using unrealistic stellar distributions and 
line emissivity. The most serious objection to the model as 
yet is the large number of supergiants needed to produce the 
line luminosity. The resulting high stellar density leads to an 
unacceptably large collisional destruction rate and to a very 
high mass injection rate into the HIM which is precluded by 
observations (Begelman & Sikura 1991). These unfavorable 
estimates depend strongly on assumptions about the stellar 
distribution in the galactic core, estimates of the BLR size 
and the size and density structure of the winds. 

In this paper we will concentrate on the questions of the 
line ratios, the number of bloated stars required for obtain- 
ing the covering factor and the size of the BLR. The issues 
relating to the dynamics of the BLR and the line profiles 
will be discussed elsewhere. In § 2 we specify the properties 
of the stellar core and of the winds of the bloated stars and 
discuss the assumptions and approximations underlying our 
models. In § 3 we describe the photoionization code, its ap- 
plication to the winds and the calculations of the emission 
line spectrum. In § 4 we discuss the various observational 
and theoretical constraints on the wind models. In § 5 we 
present results for a single bloated star and for the entire 
bloated star population. In § 6 we discuss the implications 



for the bloated stars model and the constraints that can be 
placed on it. 



2 THE STELLAR CORE AND BLOATED STAR 
MODEL 

Our aim is to model the stellar winds scenario and to calcu- 
late the expected emission line spectrum. Two ingredients 
are necessary: the stellar distribution function, n+, and the 
properties of the super giant winds. In modeling the stel- 
lar distribution function we will use the results of detailed 
numeric calculations that follow the evolution of a stellar 
core around an initial black hole and the resultant black 
hole growth and continuum radiation evolution (Murphy, 
Cohn & Durisen 1991, hereafter MCD). Since giant stars 
in AGN core may be very different from those observed in 
the solar neighborhood, we adopt the general term bloated 
stars rather than the more specific term supergiants for des- 
ignating the wind emitting objects. The bloated stars will 
be modeled by simplified wind structures with no claim to 
hydrodynamical self-consistency. The bloated stars model 
differs in several significant aspects from the clouds model: 
The wind gas density is not uniform but changes with the 
outflowing gas, The distribution of bloated stars is expected 
to be correlated to the density of the galactic core (and pos- 
sibly also to the irradiating flux) and the velocity field is 
that of the stellar core. 

2.1 The stellar core 

The attempts of previous works to use a dynamically self 
consistent stellar core distribution were based on a stellar 
distribution of the form ra* oc r -7 ^ 4 (e.g. Kazanas 1989). 
This distribution can be derived from simple dimensional 
arguments of energy conservation in the equilibrium state 
of a collisionless system of point masses in a deep potential 
well, isotropic both in position and velocity space (see e.g. 
Binney & Tremaine 1987). The immediate stellar surround- 
ing of AGN is probably isotropic and virialized irrespective 
of the host galaxy type. However, the omission of star-star 
collisions is not justified in the dense environment of the 
deep potential well. 

In this work we use the results of a much more detailed 
and realistic analysis performed by MCD. Their work follows 
the dynamical and luminosity evolution of a galactic nu- 
cleus consisting of a dense coeval stellar system surrounding 
a massive black hole. The stellar population contains stars 
of masses from 0.3 to 30Mq which are initially distributed 
according to a power-law initial mass function (IMF): 

dn^(m) oc m -1 - 1-1 "^ dm , (1) 

where dn* is the number of stars in the mass interval dm. 
The black hole is fueled by mass loss due to tidal disrup- 
tion of stars in the loss cone, inelastic (physical) collisions of 
stars, and stellar evolution. The mass loss due to collisions 
is treated in detail taking into account the radii, masses, rel- 
ative velocity and impact parameter of the colliding stars. 
However, collisions with evolved stars are not taken into ac- 
count. Mergers and the formation of binaries are also not 
considered. It is assumed that all the mass released by these 
processes accumulates in a gas reservoir near the black hole 
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and is eventually completely accreted. It is implicitly as- 
sumed that the accretion is spherical and limited by the 
Eddington accretion rate. The initial conditions are a seed 
black hole of Mbh ~ 10 4 Mq, a Plummer stellar distribution 
function with initial velocity equipartition and the power- 
law index of the IMF. The isotropy in coordinate and veloc- 
ity space is maintained throughout the calculations. 

The MCD results show that in low density cores 

7xl0 6 Mq/pc 3 ), where elastic collisions dominate, the 
core expands due to energy released by stars captured in 
tightly bound orbits and a ra* oc r -7 ^ 4 density power-law 
is established, as discussed above. In higher density cores 

7xlO r MqI pc 3 ), the inelastic collisions lead to core con- 
traction and a oc r -1 ^ 2 density power-law. In all models 
the central region is eventually depleted by stellar disrup- 
tions and collisions. This is more pronounced in high density 
cores due to the larger black hole mass and the higher fre- 
quency of collisions. At large distances (> 1 pc), where the 
stellar dynamic time-scale is very long, the density function 
essentially reflects the initial Plummer law stellar density 
and is therefore to some extent arbitrary. MCD also show 
that the density functions of the different mass groups, all 
initially the same, evolve differently ('mass segregation'). 

The calculations reproduce the full range of observed 
AGN luminosities from I0 43 to I0 48 erg/s by varying the 
initial core density and IMF power-law index. The lumi- 
nosity is calculated on the assumption of a constant mass 
conversion efficiency of O.f . Initially limited by the Edding- 
ton luminosity, it grows exponentially, reaches a sharply de- 
fined peak luminosity after a /ewxIO 8 yrs and then declines 
slowly. The black hole growth is likewise limited by the Ed- 
dington limit until the peak luminosity is reached at which 
time the accretion rate equals the replenishment rate of the 
gas reservoir. From that point on the gas reservoir is rapidly 
exhausted and the accretion of inflowing gas proceeds im- 
mediately at a rate governed by the mass loss processes in 
the core (see figures 5 & 10 in MCD). The peak luminosity 
and the final mass of the black hole increase with the initial 
core densities. 

A full study of the bloated stars scenario should in prin- 
ciple check separately each of the core configurations de- 
scribed above. Our work, being an initial study, will use one 
particular core model in the attempt to reproduce the line 
emission of a typical AGN. The choice of this core model, 
as discussed below, is an attempt to model an AGN which 
is both representative of its class and has the optimal prop- 
erties for the bloated stars scenario as suggested by the pre- 
vious works on this subject. 

Previous estimates (Scoville & Norman 1988; Kazanas 
1989; Begelman & Sikura 1991) predict that the number of 
BSs in the core may not be enough for the required BLR 
covering factor of about O.f. It is therefore preferable to 
choose a core model that has the largest number of stars 
for a given continuum luminosity. The MCD results indi- 
cate that the more massive cores have a larger fraction of 
their mass locked in the black hole and consequently fewer 
stars as compared to less massive cores. On the other hand 
less massive cores may not be luminous enough for bright 
AGNs. Our choice is to model an intermediate mass core 
that reaches bright Seyfert I or low quasar luminosities but 
still has a significant fraction of its mass in the stars. 

It is necessary to specify the core's evolutionary stage 



since the shape of the stellar distribution changes with time. 
The initial stage of exponential luminosity growth is too 
short-lived (fewxW 7 yrs) and not luminous enough. The 
choice between the two subsequent stages, peak luminosity 
and post-peak power-law decline, is less obvious. The peak 
luminosity lasts a few xfO 8 yrs at which time the stellar 
density is relatively peaked at the center. This may improve 
the covering factor problem and tend to decrease the BLR 
size (in step with recent reverberation studies) but, on the 
other hand, may also increase the collision rate. This choice 
is problematic since it implies that the lifetime of a single 
AGN is short and consequently there should be many galax- 
ies with relic supermassive black holes at their nucleus to 
account for the observed space density of AGN (Blandford 
1990). 

The alternative of modeling the core at the post-peak 
stage has the advantage that it is the stage where the AGN 
spends most of its lifetime and that the core's structure and 
luminosity change relatively slowly so that the modeling is 
not strongly time dependent^. During the post-peak stage 
the core is depleted, the collision rate is relatively low and 
the covering factor is expected to be smaller. 

In this work we chose to model a core just before its 
peak luminosity. This model (MCD model 2B) has an IMF 
index x = 1.5 (eq. I) and begins with a central density of 
7.2xIO r M / pc 3 , a core mass of 3xIO 8 M within the inner 
I pc and a total nucleus mass of 8.5xIO 8 M within 100 pc. 
The luminosity peaks at 4.0 xfO 8 yrs with L pea k = 2.4 xfO 46 
erg/s and a black hole mass of Mbh = I.9xIO 8 M . We 
model the cluster at 3x10 yrs at which time the accretion 
rate is at the Eddington limit with a bolometric luminosity 
of I0 46 erg/s, the black hole mass is 8x10 M and the mass 
loss from the core to the central gas reservoir proceeds at a 
rate of O.3M / yr. The stellar distribution is taken to be that 
of the O.8M mass bin on the assumption that stars of about 
solar mass are representative of the BSs progenitors (§ 2.2) 
(the distribution functions of the various stellar mass bins 
vary somewhat in their radial scale but their overall shapes 
are quite similar). 

The stellar density distribution exhibits the three dis- 
tinct zones discussed above: depletion near the black hole, 
a r -1 / 2 power-law at intermediate radii and steepening at 
large radii. A simplified form (Figure I) was interpolated 
between stellar distributions given by MCD at 2.8x10 and 
1.04 xfO 9 yrs and was used to approximate the stellar dis- 
tribution. The interpolation was done linearly in log(t) and 
log(ra^) since the evolution of ~ I M stars at ~ I0 8 yrs is 
governed by stellar collisions and disruption by the black 
hole, both roughly linear in the log(t) — log(M) plane. 
This distribution contains 7xIO 5 M in the inner O.fpc, 
8xIO r M in the inner I pc and a total of 2xIO 8 M which 
is ~ 1/4 of the total stellar mass of the system. 

The continuum spectrum used in this work is a syn- 
thetic spectrum typical of AGNs (e.g. Rees, Netzer & Fer- 
land 1989) with i ion = 7xf0 45 erg/s and L bo i = I0 46 erg/s 

t MCD conclude that the slow luminosity decrease of their mod- 
els is inconsistent with the observed redshift dependence of AGN 
luminosity and suggest a decreasing mass conversion efficiency as 
a way of obtaining a steeper luminosity decline. This is beyond 
the scope of the present work. 
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Distance from BH [pc] 
Figure 1. The stellar distribution function model. 

and mean ionizing photon energy of {hp) = 2.58 Ryd. We 
do not specify the origin of this continuum and assume that 
all the continuum radiation is emitted isotropically from a 
point-like source. The possibilities of an anisotropic source 
or obscuration are not considered. 

2.2 The bloated star 

BSs in the AGN may be very different from solar neighbor- 
hood supergiants. Current understanding of stellar winds 
is based on observations in the solar neighborhood and on 
theory developed for isolated spherical main sequence stars 
(Mihalas 1978). The wind mechanism of red giants and su- 
pergiants is not well understood and furthermore, there are 
yet no results about the possible effect of the AGN's in- 
tense radiation field on the wind structure. A possible mech- 
anism for driving the mass loss, internal radiation pressure 
on dust formed in the cool regions of the atmosphere (Zuck- 
erman 1980), may be strongly suppressed by the heating of 
the atmosphere by the external radiation. While the pho- 
ton flux probably can not increase significantly the mass 
loss rate (Voit & Shull 1988), a strong enough neutrino flux 
may transform low mass stars into red giants (MacDonald, 
Stanev & Biermann 1991) and the existence of a signifi- 
cant neutrino flux in AGN has been suggested (Stecker et al 
1991). 

In the absence of full reliable models, we adopt the 
known properties of SGs as the reference point. It is con- 
venient to describe the BS, like a SG, as composed of a star 
with radius and mass M + , which emits from its surface 
a wind whose denser, line emitting fraction is characterized 
by a radius R w and mass M w . 

In the presence of the unidirectional continuum radia- 
tion field, the spherical symmetry of an isolated star is re- 
placed by the cylindrical symmetry about the direction of 
the irradiating flux. For a SG of radius 10 13 cm and a 3000K 
photosphere, l/3pc from a 10 46 erg/s continuum source, the 



external flux at R w = 10 14 cm is an order of magnitude larger 
than the internal flux. The spherical symmetry must be per- 
turbed and the density and gas velocity of the wind must 
be strongly affected. In this work we will not attempt to 
deal with the anisotropy of the wind and shall assume a 
spherically symmetric wind without specifying the mecha- 
nism that is driving the wind. This can be justified if the 
results of the different wind models show that the line emis- 
sion is not very sensitive to the details of the wind's internal 
structure. 

The wind's radial profile is obtained from the continuity 
equation 

M w = ATrR 2 v(R)N(R) , (2) 

where M w is the mass loss rate, R is radial distance from the 
center of the star , v(R) is the wind velocity and N(R) is the 
gas density. In our models the wind's velocity field is just an 
alternative way of defining the density structure. Velocity 
fields which formally describe decelerating flows slower than 
the escape velocity may be better thought of as describing 
static envelopes without bulk flow. 

Voit & Shull (1988) place lower and upper limits on 
the SG fraction in the stellar population based on values 
observed in the solar neighborhood. These limits of 10 -6 
and 10 -3 respectively bracket the simple order-of-magnitude 
estimate of ~ 10 -4 derived from the mass loss rate of ~ 
lO _6 M /yr and progenitor lifetime of ~ 3xl0 9 yrs. We 
shall adopt these solar neighborhood values and assume that 
~ 0.01 of the stars in the core are red giants (RG) with 
R ~ 10i? Q and M w ~ lO" 8 M / yr and ~ 10" 4 of the stars 
are SGs with R ~ 100i? Q and M w ~ 10" 6 M Q /yr. The rest 
of the stars are main sequence stars with a typical radius 
of IRq and mass-loss rate of ~ 10 -14 MQ/yr. If we assume 
that the gas density at the boundary of the winds is the same 
for SGs and RGs (see 'Comptonization' below) and that 
the wind flows at the escape velocity v(R) = ^2GM*/R, 
then from eq. 2 it follows that the size of the wind scales as 

2 /3 

R w cx M w and the total geometric cross-section of all SGs 
is (10" 4 /10" 2 )x(10" 6 /10" 8 ) 4/3 ~ 5 times larger than that 
of the RGs. If we take the boundary density as an estimate 
for the density in the wind, the column density which scales 
as R w is ~ 20 times larger in a SG. The total mass in the 
wind scales as R% and therefore the SG winds contain ~ 100 
times more mass. These ratios justify neglecting the RG 
winds and main sequence stars and dealing only with the 
SG winds. 

The ionizing radiation may also induce mass loss from 
the star by radiation pressure at the wind base and by ab- 
lation at the star's limbs. Voit & Shull (1988) estimate this 
contribution to be negligible compared to the normal evolu- 
tionary mass loss of supergiants for all stars but those very 
near to the continuum source. In our model M w does not 
depend on the ionizing flux. We also assume that the BSs 
fraction in the stellar population is independent of the dis- 
tance from the black hole. 

The ambient HIM pressure may regulate or quench the 
stellar winds. We assume that the winds expand into vac- 
uum. Further work is needed to clarify this issue, but we 
note that Voit & Shull (1988) estimate that this assump- 
tion is valid in the somewhat different case of X-ray driven 
winds, 1 pc from the continuum source. 
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Stellar collisions are important in the dense core and are 
taken into account in modeling the core evolution (§ 2.1). It 
is reasonable to expect that the ongoing soft collisions will 
affect the wind mass and size. We, however, ignore the effects 
of these repeated collisions on the BS's structure and assume 
it remains constant throughout its life. This is equivalent to 
approximating the cumulative effect of the many collisions 
by a single destructive one. 

The radius of the line emitting inner fraction of the 
wind (henceforth the wind boundary) is an important pa- 
rameter that determines the covering factor and boundary 
(i.e. minimal) density of the wind. We consider four physical 
mechanisms that influence this radius: 

(i) Comptonization: Following Kazanas (1989), we as- 
sume that all gas with an ionization parameter larger 
than a certain critical value, corresponding to a min- 
imal boundary density Nc, goes directly to the hot 
comptonized phase and its optical depth drops to zero. 
We define the comptonization radius i?c(n>h) as the 
distance where the wind gas density, which decreases 
with R, drops to the comptonization density Nc at 
which point the gas temperature rises to ~ 5xl0 4 K. 
It should be emphasized that unlike in the 2-phase 
model, the cold and hot phases need not be in pressure 
equilibrium: the wind material is continuously lost to 
the HIM and replenished by new material from the 
star's surface. 7Vc(''bh), which is a function of con- 
tinuum flux and depends on the continuum spectrum 
and chemical composition of the gas, is calculated self- 
consistently for the model continuum(§ 3.2). 

(ii) Finite wind mass: The wind size is limited by the to- 
tal amount of mass within the wind. We define the 
average wind mass radius, i? m ass, as the radius within 
which the average wind mass is contained. The mass 
loss rate during the super-giant phase of an isolated 
star is estimated to be (Kudritzki & Reimers 1978) 

M w ~ 5.5xlO- 13 (i/i )( J R/ J R )(M /M) M Q /yr ,(3) 

with the luminosity estimated at (Paczynski 1970) 

L = 5925O(M core /M - O.478)i . (4) 

If the SG's mass is small enough it ends its life as 
a white dwarf with a mass of ~ O.6M . The mass 
loss rate is not constant during the super-giant phase 
but nevertheless we will assume that the line emitting 
gas can be represented by a wind that contains half 
of the total ejected mass. The SG mass loss rate of 
M w ~ lO -6 M /yr discussed above is consistent with 
eqs. 3 and 4 for progenitor stars of R* = 10 13 cm and 
IMq which have half of the total ejected mass, O.2M , 
in the wind and O.8M in the star's core. We use these 
values for the BSs in all the wind models presented be- 
low. A realistic population of SGs includes stars with 
both less and more than O.2M in the wind and mod- 
eling the emission lines of an average wind mass is not 
identical to modeling such cases. However, the cover- 
ing factor is the most important parameter and such 
a simplification does not alter the overall conclusions. 

(iii) Tidal disruption by the black hole: We define the tidal 
radius as i?tidai(n>h) = -X"tidai(M*/Mbh) 1/3 n>h where 
-X"tidai is a factor of order unity. Since a wind out- 



flowing at velocities above the escape velocity is not 
gravitationally bound to the star, i?tidai should be un- 
derstood in this context as the radius where the tidal 
forces disperse the gas effectively enough for immedi- 
ate comptonization. The results presented below are 
based on the assumption of Xtidai = 2. 
(iv) Maximal velocity: In some of our models, the wind 
velocity increases with R indefinitely and an artifi- 
cial maximal velocity cutoff must be introduced. This 
maximal velocity depends on the boundary conditions 
determined by the ambient interstellar gas. A value 
of 30,000 km/s was chosen on the grounds that this 
is roughly the largest velocity observed in broadened 
AGN emission lines. R v is defined as the radius where 
the wind velocity reaches the maximum velocity. In 
practice this constraint is necessary only in extreme 
cases of highly accelerated wind models. 

In all our models, the actual wind boundary radius is taken 
to be 

R w = min(i?c(rbh), i?mass, i?tidai(n>h), Rv) ■ (5) 

In certain cases when stars are very near the continuum 
source the radiation may penetrate the BS's atmosphere. 
We model the gas density of the star's atmosphere by the 
simple case of gas in thermal equilibrium in a gravitational 
potential, giving an exponential density law 

N, tra (R) = N(R*) exp (jfir(R* -R)) , R<R* (6) 

where m p is the proton mass, T ~ 3000K and g = GM*/ R* 2 
and can be taken to be constant for R ~ R+. 

Kazanas (1989) suggested that when comptonization 
determines the wind size, the boundary of the wind is not 
spherical but has a wide, comet-like tail due to the vary- 
ing angle of incidence of the irradiating flux relative to the 
radial pressure of the expanding wind. Comptonization is 
an atomic process and for a given gas chemical composition 
and continuum spectrum it is regulated only by the ratio 
between the photon and electron densities which are not 
vector quantities. For a given wind density structure, the 
comptonized boundary in the wind hemisphere that faces 
the continuum is therefore simply a surface of constant den- 
sity. The boundary surface of the back hemisphere depends 
on the front-to-back optical depth. In the limbs, where it is 
negligible, the boundary still follows a constant density sur- 
face. At the center, the exact shape of the wind is irrelevant 
to the line emission as the optical depth is infinite. Between 
these two limits the back hemisphere's shape may deviate 
from the surface of constant density. In this work we ignore 
this possibility and conclude, unlike Kazanas (1989), that 
a spherically symmetric density structure implies a spher- 
ical comptonized boundary (it is possible that a realistic 
wind density model, such that takes into account the mo- 
mentum and heat imparted by the continuum, will indeed 
have an elongated or comet-like comptonized boundary, but 
this would be a consequence of the non-spherical surfaces of 
constant density, and not of the flux' angle of incidence). 

In the modeling we ignore the possible obscuration of 
one BS by another. This is consistent with the assumption 
that the BLR covering is no larger than ~ 0.1. 
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3 CALCULATIONS 

3.1 Photoionization of stellar winds 

The photoionization code used in this work is ION (Rees, 
Netzer & Ferland 1989; Netzer 1993) which implements an 
iterative algorithm that is based on the approach of treat- 
ing the atomic physics in detail while approximating the 
radiative transfer. The treatment of the radiative transfer 
assumes that the gas has an infinite slab geometry which is 
homogeneous in density in the transverse plane but may vary 
in the longitudinal direction. ION receives as input the gas 
chemical composition, the longitudinal density profile, the 
external radiation spectral shape and luminosity and pro- 
duces as output the emission line luminosities. Only "solar 
composition" models, given by 

(He : C : N : O : Ne : Mg : Si : S : Fe)/H = O.OOOlx (7) 
(1000 : 3.7 : 1.1 : 8.0 : 1.1 : 0.37 : 0.35 : 0.16 : 0.4) 

are used in this work (see § 5.1.4). 

The application of ION to the stellar winds poses some 
difficulties. The geometry of the wind is far from slab-like. 
As will be described in § 3.2, the actual calculations assume 
a cylindrical shell geometry with a density gradient in both 
the transverse and longitudinal directions. Another problem 
particular to stellar winds is the existence of a velocity gra- 
dient due to acceleration. For high enough accelerations or 
lines of relatively small optical depths, the Doppler shift of 
the line profiles due to the bulk motion of the wind mass 
may affect resonance scattering and fluorescence. While this 
can be handled by ION, it was not implemented for the sake 
of simplicity and because the uncertainties due to the wind 
structure dominate the results. Various tests that were car- 
ried out suggest that two types of lines are most affected. 
The intensity of the high excitation, temperature sensitive 
lines, such as OviA1035, decreases with the increase in line 
width because the reduced optical depth in La and other 
major coolants causes a drop in the electron temperature. 
Intercombination lines, especially ClIl]A1909, can increase 
their intensity due to the fact that their optical depth is, in 
some cases, non-negligible. The increased line width reduces 
the line optical depth thus raising the critical de-excitation 
density. Both these types of lines are not affected by a large 
factor. 



3.2 The numeric procedure 

Calculating the line emission from a given core and wind 
model requires that the line emitting gas be divided into 
roughly slab-like sections such that can be handled by ION. 
This proceeds as follows: The stellar distribution is binned 
into concentric spherical shells. The BSs within a given shell 
are represented by a single BS whose wind size and density 
are typical of the shell's mean distance from the black hole. 
The wind is then divided into concentric annular cylindrical 
slices whose axis is parallel to the flux of the irradiating 
continuum. The line emission from each slice is calculated 
separately by ION given the depth dependent density profile 
and the results are then integrated back to give the total 
contribution from the wind of a single BS in the shell. The 
total line emission from the core is calculated by adding the 



contributions from each shell, weighted by the number of 
BSs in it. 

The stellar core is cut into a set of radii {r t } between 
r; n and r ou t such that: 



r,+i 



and 



^ ^ratio 



r n+ dr 

J I in 



r 2 ra+ dr 



(8) 



(9) 



where r ra tio and ra ra tio are parameters chosen to ensure that 
the core is finely binned. The values r ra tio = 1-25 and 
Wratio = 1.0 (i.e. no constraint on the fraction of stars in 
the shell) were found to give satisfactory results. 



Our photoionization calculations show that r; n 



10 1 



cm is the nearest distance where a sharp transition zone 
can still be established between the comptonized and cold 
phases as the high gas density leads to collisional suppression 
of the coolants. r ou t is chosen so as to be large enough for 
the core's integrated line ratios to approach their asymptotic 



values. The numerical results indicate that r ou t 



is a suitable value. The core is typically divided into ~ 45 
spherical bins. 

The wind's density structure is parameterized by the 
functional form of v(R), which is assumed to be indepen- 
dent of the distance from the black hole, and by the wind 
size R w . The calculation of R w (r) from eq. 5 requires the 
comptonization radius Rc and therefore the boundary den- 
sity Nc- Nc(r) is calculated self-consistently and interpo- 
lated from a table prepared beforehand for the specific con- 
tinuum spectrum used here (§ 5.1.2). The table is prepared 
for a range of rbh values by applying ION to a gas slab of 
constant density and by varying this density until the cold 
phase is established very near the slab's illuminated surface. 

The binning of the irradiated cross-section of the wind 
proceeds by cutting it into a succession of annular slices 
with radii {R t } from R = to R w such that the slice's area 
relative to the wind's cross-section satisfies 



p2 p2 

R 2 



< Rratio 



(10) 



and the maximal wind density along the i'th cut, N t , and 
that of the i + l'th cut, 7V 8 +i satisfy 



(11) 



where y4. ra tio and 7V ra tio are given parameters. The first con- 
straint is intended to ensure that the binning of the cross- 
section is fine enough. The second is intended to minimize 
the transverse density gradient in the gas slices. The values 
Rratio = 0.15 and ra ra tio = 0.7 were adopted by trial-and- 
error. This choice, which results in very narrow slices (width 
to depth ratio < 0.001 in certain cases) in contradiction 
with the infinite slab assumption, represents a compromise 
between conflicting requirements. The number of slices per 
star depends on the density gradient in the wind and ranges 
from 7 to ~ 60. The luminosity in a given line I from a single 
BS, Lt(ri), is integrated by summing the line flux from each 
slice weighted by the slice's cross-section. Lt is interpolated 
to any rbh by cubic spline. 
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Table 1. A composite quasar and Seyfert 1 observed line spec- 
trum. 



Line type Ft/F^a range F{/Fi, a range 

Line BLR only BLR and NLR 



H and He 



Ha A6563 


0.5 


± 


0. 


.3 


0.5 


± 


0. 


.3 


H/9 A4861 


0.115 


zc 


0. 


.085 


0.115 


4- 
zc 


0. 


.085 


HelA5876 


0.0175 


± 


0. 


.0125 


0.0175 


zb 


0. 


.0125 


HellA1640 


0.09 


± 


0. 


.06 


0.09 


zb 


0. 


.06 


Resonance 


















Civ A1549 


0.55 


± 


0. 


.25 


0.55 


zb 


0. 


.25 


IN V AliiU 


0.225 


4- 
zc 


0. 


.175 


0.225 


4- 
zc 


0. 


.175 


o vt ). i o 3 5 


0.375 


4- 
zc 


0. 


.325 


0.375 


4- 

zc 


0. 


.325 


\/ f^v\r hi ff h PYntatinn 
VC1V llllill CA^lldllUll 


















CiiiA977 


0.05 


± 


0. 


.05 


0.05 


zb 


0. 


.05 


Intercombination 


















Cm] A1909 


0.275 


± 


0. 


.125 


0.275 


zb 


0. 


.125 


Low excitation 


















MgllA2798 


0.275 


± 


0. 


.225 


0.275 


zb 


0. 


.225 


Forbidden 


















[O ill] A4363 


0.005 


± 


0. 


.005 


0.015 


zb 


0. 


.015 


[O ill] A5007 


0.0075 


± 


0. 


.0075 


0.5 


zb 


0. 


.5 


[Nev] A3426 


0.005 


± 


0. 


.005 


0.225 


zb 


0. 


.175 


[Fexi] A7892 


0.005 


± 


0. 


.005 


0.01 


zb 


0. 


.01 



4 OBSERVATIONAL AND THEORETICAL 
CONSTRAINTS 

4.1 The observed emission line spectrum 

We calculate the emission in 15 lines (table 1) as well as 
the total reprocessed emission in lines and continuum for 
each of our models. These lines were chosen for the rel- 
ative reliability in their calculation, the statistical signifi- 
cance of the observed data and their value as probes of the 
conditions in the gas. The hydrogen and helium lines La, 
Ha, H/9, HelA5876 and Hell A1640 are useful for estimating 
the total absorbed flux and the covering factor. The reso- 
nance lines OviA1035, NvA1240 and CivA1549 are good 
indicators of the gas temperature and ionization parame- 
ter and the intercombination line C III] A1909 is an indicator 
of intermediate densities of ~ 10 9 cm -3 . The low excita- 
tion line Mg II A2798 is an indicator of cold, partly ionized 
gas. The very high excitation line C III A977 is an indicator 
of high density, high temperature gas. The forbidden lines 
[Nev] A3426, [Ora]A4363, [Ora] A5007 and [Fexi]A7892 are 
indicators of low density gas at various levels of excitation. 
The Ov]A1218 line underlying the La may introduce, in 
principal, an uncertainty in the observed La flux. However, 
since our calculations show that the typical Ov]A1218/La 
ratio in our wind models is less than 1%, we omit this line 
from the list above. 

The calculated line emission from the wind models is 
expressed in terms of the La emission and compared to a 



composite quasar and Seyfert 1 emission line spectrum given 
in table 1. The range around the mean observed value was 
chosen so as to contain about 95% of the observed objects. 
These values are intended for reference only and do not rep- 
resent a sample mean and standard deviation in the rigorous 
statistical sense: For some lines the scatter around the mean 
is actually not symmetric and the 95% range implies in fact 
a zb2<r rather than zblcr interval. A convenient standard for 
expressing the discrepancy between a calculated line I and 
the observed range is in terms of the normalized difference, 
Li = (Lt/L-Lc — fit) /at, where fit and at are taken to be 
the mean and scatter in table 1. In addition to the diag- 
nostics provided by the individual lines, a \ 2 score is calcu- 
lated with these fit and at values as a measure of the overall 
goodness-of-fit of all the calculated line ratios. The number 
of degrees of freedom is 12 = 14 — 2, the number of indepen- 
dent line ratios minus the number of free parameters in the 
wind density profile. Since Lt formally resembles the stan- 
dard deviation difference, we shall use below terms such as 
'3cr discrepancy', but it should be emphasized that in view 
of the approximate nature of fit and at, the calculated val- 
ues of Li and \ 2 should be interpreted only as rough relative 
indicators of the models' success. 

Reddening is a possible complication in comparing the 
calculated line ratios to the observed ones. The issue of the 
intrinsic reddening in AGN is still unclear (MacAlpine 1985). 
The motivation for assuming some intrinsic reddening is the 
very large deviation between the observed and calculated 
hydrogen line ratios, as well as the ratios of weaker lines 
such as O I and Hell (Netzer 1990). In the framework of the 
screen hypothesis, the application of reddening corrections 
to these deviations typically gives -Eb-v values of 0.05 to 
0.25 with a mean of 0.15. Here we use an analytic approx- 
imation of the extinction curve based on interstellar obser- 
vations (Seaton 1979). We correct the results of each wind 
model, integrated over all the BSs in the core, so as to ob- 
tain the mean observed La/H/5 1 ratio and then correct all the 
other lines accordingly. In this scheme, similar to that used 
to deduce the intrinsic reddening in cloud models, -Eb-v is 
not an additional free parameter and may have a different 
value for each of the wind models. We do not consider the 
possibility of reddening by internal dust. 

4.2 The covering factor and fraction of BSs 

The BLR covering factor, Cf, can be translated into the 
number of BSs once their line emission is calculated. A pos- 
sible definition for the covering factor in terms of the AGN's 
spectral properties is the ratio between the total reprocessed 
BLR luminosity to the total continuum luminosity, in which 
case the canonical value is Cf ~ 0.1. Alternatively it can 
be defined in terms of the observed La equivalent width in 
which case the canonical value is _EWBLR(La) ~ 100 A. We 
choose the latter since it is more closely related to the quan- 
tities directly observed. The required fraction of the of BSs 
in the stellar population, /bs , is given by 



J^ ut 4:7rr 2 n4r)EW^La)(r)dr ' 

where EW^(La) is the contribution of a single BS to the 
La equivalent width and we assume that the fraction of the 
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BSs does not depend on rbh. The choice of the integration 
limits r; n and r ou t will be discussed below in § 4.4. Since 
ji* represents only the low mass stars which account for 
~ 1/4 of the stars in the core (§ 2.1), /bs can take in the 
extreme case (all the stars are bloated) the maximal value 
~ 4. If solar neighborhood values apply we expect /bs ~ 
10 -4 . Wind models with /bs ^ 1 imply inconsistency with 
the assumed stellar core model. 

4.3 Mass loss constraints 

An important constraint on the number of BSs is mass 
loss due to stellar evolution in the bloated phase and col- 
lisions between BSs. These processes are not part of the 
MCD calculations which do not assume the existence of BSs. 
Mass loss from main sequence stars is already included in 
the 0.3MQ/yr rate calculated by MCD for our core model 
(§ 2.1). BS mass loss significantly in excess of this value is 
inconsistent with the core model and may also imply a large 
optical depth to electron scattering, r es > 1, contrary to ob- 
servations. We will use here a simple model for the mass loss 
due to collisions between two BSs: The average fractional 
mass lost at each collision, k, is about 0.1 for a collision of 
two equal mass stars approaching each other with infinite 
velocity. The fractional mass loss decreases by up to two 
orders of magnitude for slower collisions and for stars with 
different masses (Murphy, Cohn & Durisen 1991). We there- 
fore ignore MS-BS collisions and assume k = 0.1 for BS-BS 
and MS-MS collisions. We estimate the velocity dispersion 
in the core by 



<r*(r) = ^/GM(r)lr 



(13) 



where M(r) is the total mass (black hole and stars) within 
radius r. The collisional mass loss rate per volume due to 
BS-BS collisions is 

PBScoii(r) = (/Bsra*(r)) 2 (T^(r)K7rx (14) 
[(M* - M w )Rl + M w (R 2 w (r) - , 

where the two different cases of collisions between the BSs 
themselves and between the winds are considered separately. 
We also estimate the mass loss due to MS-MS collisions, 
independently of the MCD result of 0.3MQ/yr, by assuming 
that the the 1 to 4 overall ratio between and the total 
MS stars in the core holds at every rbh- Using solar values 
to represent the MS stars we obtain 

PMScoii(r) = (max(0, 4 - /Bs)ra*(r)) 2 <T*(r)K7rM _R 2 . (15) 
The evolutionary mass loss rate per volume from BSs is 



Pevol = Af w / B sw*( r ) 



(16) 



These estimates implicitly assume that the ejected mass is 
spread evenly over the spherical shell containing the star's 
orbit. The total mass loss rate from these two processes is 
/•rout 2 

M B s = / 47rr (pBScoii + pMScoii + pevoi) dr . (17) 

J r\ n 

On the assumption that this mass is neither accreted nor 
ejected from the AGN, the upper limit on the buildup rate 
of the optical depth to electron scattering is 

Tes = <JT (pBScoll + pMScoll + Pevol) dr , (18) 

J Fin 



where <tt is the Thompson cross-section. Since the BLR size 
is of the order of at least 0.1 light year, instant removal of 
this mass is not possible. Even if all this material is even- 
tually accreted, f es xlyr provides an estimate for the mean 
electron scattering opacity during the phase of the black hole 
buildup. 

4.4 BLR size constraints 

The integrated luminosities, the number of BSs, the BS frac- 
tion in the stellar population, the magnitude of the mass loss 
effects and the size of the line emitting region all depend on 
r; n and r ou t and therefore constrain their values. The ab- 
sence of an intermediate zone between the BLR and NLR 
may mean that there is some upper cutoff mechanism to the 
BLR, e.g. dust (Netzer & Laor 1993). While it may be pos- 
sible that the decreasing line emission efficiency of the BSs 
will provide a natural explanation for this cutoff, it is also 
possible that our simplistic assumptions that M w and the 
fraction of BSs in the population are independent of rbh are 
wrong and an artificial cutoff must be introduced to correct 
for it. In this work we calculate the line emission up to a 
distance of 56 pc (§ 3.2), which is well beyond the estimated 
BLR size, in order to check this matter explicitly. The exact 
choice of r; n should not affect the the overall results strongly 
as there are very few stars at the inner radii (although in cer- 
tain cases their contribution may be significant — see § 5.2.2). 
In this work we will assume that r; n = 2 x 10 16 cm (§ 3.2) and 
treat r ou t as yet another free parameter of the model. 

As is apparent from eq. 12, /bs is a decreasing function 
of r<,ut. The number of BSs participating in the BLR line 
emission, Nbs, is given by 



^VBs(r ou t) 



100 A f r ° ut r 2 njr)dr 

J 'in 

/ r j° u ' r 2 n4r)EW^La)(r)dr 



(19) 



and is not necessarily a rising function of r ou t . 

Observational information on the BLR size obtained 
from line reverberation campaigns gives a rough estimate 
of 

1/2 



rBLR 



0.1 



L 



10 46 erg/s 



pc 



(20) 



Here we calculate only the flux weighted radius of a line 
irir' r3n *( r ) L e(r) dr 



J r; ° ut r 2 n*(r)L e (r) dr 



(21) 



We will not perform a detailed comparison between the flux 
weighted radius of different lines and limit ourselves to rough 
consistency checks between the rf w of the prominent BLR 
lines and rBLR of eq. 20. 



5 RESULTS 

The viability of the bloated stars scenario depends critically 
on the properties of the single BS as a line emitting object. 
In presenting the results we aim to separate issues that can 
be analyzed by studying the properties of the single bloated 
star from those that depend on the properties of the entire 
bloated star population in the AGN core. Following that, we 
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try to isolate those features of the bloated stars that do not 
depend strongly on the choice of wind structure and explore 
to what extent the structure dependent properties can be 
varied. 



5.1 The single bloated star 

5.1.1 The wind models 

The wind model that was considered in previous works 
(Kazanas 1989; Begelman & Sikura 1991) was that of a con- 
stant velocity flow with v ~ 10 6 cm/s. Here we explore a 
much wider range of wind structures with two forms of ve- 
locity profiles 

(i) A power-law profile 

v(R) = v^R/R*)-" (22) 

where v+ is the velocity at the base of the wind at R*. 
This form includes the case of constant velocity, a = 
0, and of free-fall motion, a = 0.5. We explore this 
form in the range a = —2.0 to 0.5 and v* = 5xl0 3 to 
5xl0 6 cm/s. The constant velocity case (a = 0) with 
v* = 10 6 cm/s will serve below as the reference model. 

(ii) An asymptotic profile 

v(R) = v 00 (l-(R/R*)- p ) + v* (23) 

where Vco is the velocity at infinity. This analytic ap- 
proximation reproduces the main features of Voit and 
S hull's (1988) numeric hydrodynamical work on radi- 
ation driven winds, v* = 5xl0 4 cm/s is the critical 
velocity at the base of the wind which is assumed to be 
a constant independent of Vco and fi. We explored this 
form in the range fi = 0.5 to 1.0 and «oo = 5xl0 5 to 
5xl0 6 cm/s. 

Both velocity profiles have two free parameters: a and v+ 
for power-law winds and fi and «oo for asymptotic winds. 

The asymptotic flow approaches constant velocity at 
large radii and is therefore expected to produce an emission 
line spectrum similar to that of a constant velocity wind 
of the same velocity. Figure 2 presents a comparison of the 
emission in several representative lines from a constant ve- 
locity star of v = 10 6 cm/s and and an asymptotic wind flow 
(eq. 23) with «oo = 10 6 and fi = 0.75. The boundary density 
of both wind models, Nc, is determined by comptonization 
and is therefore the same for both at all distances. As ex- 
pected, at rbh > 10 17 cm, where i? w ^> R+, the differences 
in the line emission relative to La is no more than 10% 
and are well within the numeric uncertainty. At smaller dis- 
tances the constant wind model exhibits emission typical of 
colder, denser gas relative to the that of the asymptotic wind 
model: weaker forbidden and semi-forbidden lines, lower 
CmA977and Cm] A1909emission and higher Nv and Mg II 
emission. 

These results can be explained by the difference between 
the density gradient in the two models: The asymptotic wind 
gas is denser near (a consequence of its lower wind base 
velocity v* = 5xl0 4 cm/s) and therefore it is comptonized 
at R w > whereas the lower gas density of the constant 
velocity wind is comptonized at R w < R+. The density gra- 
dient of the constant velocity model is the steep exponential 



r N MgII 2798 

' \ 
' \ 
' \ 




.18 



r BH [cm] 

Figure 2. Ratios of line emission in two wind models: a constant 
velocity wind of = 10 6 cm/s and an asymptotic now wind of 
v(K) = 5X10 4 + 10 6 (1 - (i?/10 13 )" 75 )cm/s. The ratio plotted 
is between unreddened luminosities in units of the ha line. 



Table 2. Boundary density, temperature and ionization param- 
eter as functions of distance from the continuum source. 



log(r 6 ft/cra) 


log(7V c /cm- 3 ) 


T C K 


U 


16.0 


12.1 


7xl0 4 


2.8 


16.5 


11.0 


7xl0 4 


3.5 


17.0 


10.0 


7xl0 4 


3.5 


17.5 


9.0 


7xl0 4 


3.5 


18.0 


8.0 


7xl0 4 


3.5 


18.5 


7.0 


8xl0 4 


3.5 


19.0 


6.1 


9xl0 4 


2.8 


19.5 


5.2 


9xl0 4 


2.2 


20.0 


4.2 


5xl0 4 


2.2 


20.5 


3.2 


5xl0 4 


2.2 


21.0 


2.2 


5xl0 4 


2.2 


21.5 


1.2 


5xl0 4 


2.2 



gradient of the stellar atmosphere which does not allow the 
existence of an extensive hot, highly ionized zone. 

Similar results were also obtained for asymptotic wind 
models with fi = 0.5 and 1.0 and Voo ranging from 5xl0 5 to 
5xl0 6 cm/s. We conclude that the asymptotic wind models 
do not differ significantly from constant velocity models in 
their line emission spectrum and in the following we will 
concentrate our study on the power-law velocity models. 

5.1.2 Comptonization and limits on broad forbidden lines 

Under the assumption of comptonization, the outer bound- 
ary of the line emitting region of the wind is defined by the 
onset of instabilities in the cold gas which drive it to the 
hot phase where it becomes part of the interstellar medium 
(Kazanas 1989). Unlike previous works, we calculated the 
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instability conditions self-consistently for our assumed con- 
tinuum. We find that the transition to the unstable branch 
occurs at a temperature of about 50,000K. The line emis- 
sion properties of the gas at the comptonization radius are 
determined by the ionization parameter, 



U 



f °° L v Ihv dv 
4irr 2 cN e 



(24) 



Table 2 displays the dependence of the boundary density, 
temperature and ionization parameter on the distance from 
the continuum source. We note that unlike the 2-phase 
model which predicts that Nc oc f/r^h an d that U is 
constant, our calculations show that Nc actually decreases 
somewhat slower than f/r^h an d U decreases with the dis- 
tance. Of greater significance is the fact the values of U at 
the onset of the stable cold phase are 4-6 times higher than 
the BLR range of U = 0.3 to f inferred from photoioniza- 
tion studies and the one implied by the boundary densities 
assumed by Kazanas (4989). 

Gas with ionization parameter values as high as those 
of our model has the distinct signature of strong very high 
excitation lines. In particular, we find that some high ex- 
citation forbidden lines, with large critical densities (e.g. 
[FeXl] A7892) can dominate the emission in the outer parts 
of some winds. Although the wind density is not uniform 
but increases inwards, the spherical wind is dominated by 
the thinner limbs. This is a simple consequence of the wind 
geometry: An infinite slab of thickness R and a thin outer 
layer of thickness R — R' , such as commonly used in numeri- 
cal modeling of BLR clouds, has only 4 — R' / R of its volume 
occupied by the low density layer, whereas in a sphere of ra- 
dius R the same layer occupies 4 — (i?'/i?) 3 of the volume. 
As a consequence, the total emission line spectrum of the 
wind is strongly biased towards the high ionization parame- 
ter spectrum (figure 3). Figure 4A shows that the forbidden 
lines emission from the reference wind model far exceeds the 
observed BLR values as long as the density at the comp- 
tonization radius is lower than the critical density of the 
forbidden line. In particular, the [FeXl]A7892 emission re- 
mains high (> +3<t) until as close as ~ 0.2 pc. Similar results 
were obtained for all wind models whenever the process that 
determines the wind boundary is comptonization. 

We considered two alternatives for lowering the ioniza- 
tion parameter and thereby for relaxing the constraints on 
the comptonization mechanism. The first is the assumption 
of a softer continuum spectrum with mean ionizing pho- 
ton energy of {hp) = 4.27 Ryd. This spectrum lowers U 
to 2.4, a value still too high for the BLR, and therefore 
did not change the forbidden lines luminosities significantly 
(figure 4B). The second alternative is the assumption of a 
lower critical temperature of ~ 3x40 4 K. A possible justifi- 
cation may be that the instabilities between the two phases 
preclude an equilibrium very near the cross-over tempera- 
ture. This critical temperature lowers U to 0.7-0.9, thereby 
increasing the collisional suppression distance of the forbid- 
den lines to ~ 4 pc, at which point the [Ne v] A3426 emission 
increases above +3<r (figure 4C). 



K 

V 



y-^FeXl] 7892 




[OIII] 5007 



0.2 



0.4 



0.6 



0.8 



R/Rw 



Figure 3. Cumulative forbidden line emission from the wind's 
irradiated cross-section for a star at r^h = 1 pc and a constant 
velocity wind with v = 10 6 cm/s. The emission is given in stan- 
dard deviations from the BLR values of table 1, thus all numbers 
greater than ~ 2 imply excess emission. R/R w = is pole-on and 
1 is at the outer boundary. 



5.1.3 Dependence on wind structure 

Table 3 compares the line emission from power-law wind 
models of different indices and wind base velocities at 0.4 pc 
and 4pc. The wind boundary in all models is determined 
by comptonization and therefore the boundary density, Nc, 
is identical in all models at a given distance from the black 
hole. N c = 40 9 cm" 3 at 0.4 pc and N c = 40 r cm" 3 at 4pc. 
The dependence of size of the wind on a and v+ is obtained 
from the continuity equation 2: 



R* 



M 



AttNcv^ 



l/(2-c) 



R, 



\Nc) 



(25) 



where N* is the gas density at the base of the wind. 

The constant velocity models (a = 0) have wind base 
velocities and densities spanning 3 decades, and therefore 
wind sizes and mean density gradients spanning 4.5 decades. 
The forbidden lines [O III] A4363 and [O III] A5007 are the 
most sensitive to variations of v+, decreasing by more than 
a factor of 8 at 0.4 pc and by factors of 5 and 7 respec- 
tively at 4 pc as v* increases from 5x40 3 to 5x40 6 cm/s. 
The semi-forbidden line C III] A4909 decreases with velocity 
by a factor of 5 at 0.4 pc but by less than a factor of 2 at 
4 pc. The Mg II A2798 line, much lower than its observed 
value (§ 5.4.4), increases with the velocity by up to a factor 
of 4 at 0.4 pc but is almost constant at 4 pc. The N V A4240 
line increases with the velocity by a factor more than 2 at 
0.4 pc and by a factor of 3 at 4pc. The ClIlA977 decreases 
with velocity by more than a factor of 2 at 0.4 and 4 pc. 
All the other lines change by less than a factor of 2. Some 
general trends in the results can be readily understood in 
terms of the changing ratio between the high density and 
low density gas volumes: The decrease of R w with v+ lowers 
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Figure 4. Total forbidden line emission from the wind of a single 
star, L, as function of distance from the black hole for a constant 
velocity wind v = 10 6 cm/s. The emission is given in standard 
deviations from the BLR values in table 1. A) Standard ionizing 
continuum {{hp) = 7.74 Ryd) with a critical boundary tempera- 
ture of ~ 6xl0 4 K. B) Soft continuum of (hv) = 1.27Ryd with 
critical boundary temperature of ~ 6xl0 4 . C) Standard ionizing 
continuum {{hp) = 7.74 Ryd) with a critical boundary tempera- 
ture of 3X10 4 K. 



the fraction of high density gas in the bloated star with the 
consequence of decreasing the forbidden and semi-forbidden 
lines and increasing the Mg II A2798 line. The greater over- 
all sensitivity of the line spectrum to changes in the velocity 
at 0.1 pc relative to that at 1 pc has a similar explanation: 
Although R w decreases with the velocity by 10 1 5 at both 
0.1 and 1 pc, it is ~ 10 times larger at 1 pc whereas 
is constant. The fraction of the high density gas at 1 pc is 
therefore much smaller and the wind is dominated by the 
thin outer layers regardless of v+. 

The accelerated flow models all have the same base ve- 
locity of v* = 10 6 cm/s and base density = 3xl0 10 cm -3 . 
The dependence of the size of the wind on Ac (eq. 25) is 
such that R w increases with the wind's power-law index a 
by a factor of 25 at 1 pc and only by a factor of 4 at 0.1 pc 
as a increases from —2.0 to 0.5. Similarly, vc, the velocity 
at the comptonized boundary, decreases with a by a factor 
of 670 at 1 pc and only by a factor of 16 at 0.1 pc. The 
mean density gradient in the wind changes by a factor of 7 
at 0.1 pc and by a factor of 30 at 1 pc. The variations in 
the wind structure with a are smaller than that of the con- 
stant velocity models, and this is reflected in the behavior 
of the emission line spectrum. Although the line spectrum 
changes with decreasing wind radius similarly to that of the 
constant velocity models, the magnitude of the changes is 
smaller and is more pronounced at 1 pc rather than at 0.1 pc 
due to the much greater variation range. The forbidden lines 
[O III] A4363 and [O III] A5007 increase with a by a factor of 
2.5 at 0.1 pc and by factors of 4 and 6 respectively at 1 pc. 
The forbidden [Nev]A3426 line increases with a by more 
than a factor of 2 at 0.1 pc and by a factor of 4 at 1 pc. The 
forbidden [FeXl] A7892 which changes by less than a factor 
of 2 at 0.1 pc decreases with a by a factor of 4 at 1 pc. The 
Mg II A2798 line increases with a by a factor of more than 2 
at 0.1 pc and 3 at 1 pc. The OviA1035 line is almost inde- 
pendent of a at 0.1 pc but increases with a by a factor of 3 
at 1 pc. All other lines change by less than a factor of 2. 



5.1.4 The MgllA2798/La and N v A1240/La deficiency 

The two line ratios MgIlA2798/La and NvA1240/La are 
too low compared to the observed BLR values regardless of 
the wind model parameters or the distance from the black 
hole. The Mg II A2798/La deficiency is the most conspicuous 
discrepancy. This discrepancy is exhibited by all the power- 
law wind models at distances less than ~ 1 pc from the black 
hole and up to distances of ~ 20 pc by all winds but those of 
dense and slow decelerating flows. Reddening can improve 
the ratio only by a factor of ~ 2 even for the extreme case of 
Eb—v = 0-2. Figure 5 shows the the range of the Mgll/La 
ratios emitted by a single bloated star from power-law wind 
models covering the full parameter range as function of the 
distance from the black hole. The Mg Il/La ratio is less than 
0.03 {L = —1.1) for all wind models up to ~ 1 pc and less 
than 0.07 {L = —0.9) up to 56 pc for all but the densest and 
slowest wind models. 

The Mg II is an indicator of cold, partly neutral gas. The 
Mg II emission pattern reflects an interplay between the wind 
temperature and column density: At small distances from 
the black hole, a significant fraction of the wind cross-section 
is occupied by the bloated star itself which can provide the 
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Table 3. The dependence of the line emission spectrum of a single BS on the power-law wind parameters. 
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Top left panel: constant velocity winds with different wind base velocities at 0.1 pc. Top right panel: as in top left panel, at 1 pc. 
Bottom left panel: accelerated flow winds with the same wind base velocity at 0.1 pc. Bottom right panel: as in bottom left panel, at 1 
pc. The line luminosities are given relative to La, unreddened. 

high column densities required for cold, partly ionized gas. 
However, the low ionization is attained at high column den- 
sities where the optical depth of Mg II is very large. At larger 
distances, where the wind is dominated by the low density 
outflowing gas and r is low, the relatively high ionization pa- 
rameter maintains the gas at a state where the magnesium 
is ionized several times. The combination of low ionization 
and low r is found only in the very slow and dense wind 
models, whose size, and therefore boundary density, are de- 
termined by the mass constraint i? m ass (Fig. 6) rather than 
by comptonization. This results in a much higher boundary 
density and consequently a much lower ionization parameter 
and low gas temperature. 

We investigated whether it is possible to increase the 
Mg II emission in the comptonized reference wind model by 



using the softer continuum or by lowering the boundary tem- 
perature as in § 5.1.2. this resulted in an increase of less than 
a factor of 2 in the former and less than a factor of 3 in the 
latter and therefore does not change the general conclusion 
that the comptonized wind models show a Mg II deficiency. 

The deficiency in the N V A1240 /La ratio is another per- 
sistent feature of the results. Again, this deficiency cannot 
be explained by reddening. Adopting a lower comptoniza- 
tion temperature increased the ratio only at rbh < 0.1 pc by 
less than 20%. The softer continuum luminosity increased 
this ratio by only ~ 20% at all radii which is still very low 
compared to the observed values. 

The NvA1240/La problem is related to the general 
composition problem discussed recently by Hamann & Fer- 
land (1993). These authors suggest that the abundance of 
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Figure 5. The Mgll(rbh) / L^(^bh) luminosity ratio from power- 
law wind models. Solid line: a = 0.5, dotted line: a = 0.0, short 
dashed line: a = —1.0, long dashed line: a = —2.0. 
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Figure 7. EW(La) per bloated star as function of the distance 
from the black hole. Solid line: a = 0.5, dotted line: a = 0.0, 
short dashed line: a = —1.0, long dashed line: a = —2.0. 




Figure 6. Wind sizes for power-law wind models. The wind 
boundary is determined by comptonization subject to the con- 
straints of X t idal = 2, Mw < 0.2Mq and v mSLX = 0.1 c. Solid line: 
a = 0.5, dotted line: a = 0.0, short dashed line: a = —1.0, long 
dashed line: a = —2.0. 



nitrogen and other elements in high luminosity, high redshift 
quasars is much larger than their solar abundances, causing 
a significant increase in the above ratio. If this is indeed 
the case, our solar composition models should not be com- 
pared to the average NvA1240/La ratio, which is affected 
by the inclusion of such quasars in the sample. We did not 
investigate high metalicity models and leave this possible 
discrepancy for future work. 



5.1.5 Line emission efficiency of bloated stars 

The line emission efficiency of the bloated star depends on 
both its size and density structure. Figure 6 shows the wind 
sizes of the power-law models as function of the distance 
from the black hole. The wind size is calculated from eq. 5. 
The tidal disruption constraint affects only the slowest and 
densest wind model of a = 0.5, v* = 5xl0 3 cm/s and the 
maximal velocity constraint only the fastest and most accel- 
erated wind model of a = —2.0 and v* = 5xl0 6 . The size of 
the slow wind models is determined at the larger distances 
by the constraint on the mass in the wind and the denser 
the wind the earlier the onset of the wind mass constraint. 
As rbh decreases, the sizes of the winds decrease until they 
reach R* = 10 13 cm, the size of the bloated star itself. From 
that point on the wind size remains nearly constant due to 
the very steep exponential density profile of the stellar at- 
mosphere. The distance at which R w = decreases as the 
wind density increases. 

Figure 7 shows the contribution of a single star to the 
La equivalent width as function of the distance from the 
black hole. This contribution is less than 10 -4 A at all dis- 
tances for all but the slow and dense models. The peaks of 
the contribution coincide with the distance where the wind 
size stops increasing with Rc and is constrained by i? m ass • 
From that point on the irradiating flux decreases as l/r^ h 
and the contribution drops steeply. 

Previous estimates of the wind's actual line emission 
efficiency were based on the assumption that the wind's ge- 
ometric cross-section, or a constant fraction thereof, is op- 
tically thick. Figure 8A shows the ratio between the total 
reprocessed luminosity from the wind to the total ionizing 
luminosity that is irradiating the wind's geometrical cross- 
section wRyj. At small distances the ratio exceeds 1 due to 
the large column densities of the bloated star itself which 
also absorb photons redwards of the Lyman limit. As rbh 
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Figure 8. Effective covering factor relative to the geometrical 
cross-section of the power-law wind models as function of the dis- 
tance from the black hole. A) Expressed in the ratio of reprocessed 
to irradiating luminosity. B) Expressed in the ratio of La photons 
to infalling irradiating ionizing photons. Solid line: a = 0.5, dot- 
ted line: a = 0.0, short dashed line: a = —1.0, long dashed line: 
a = -2.0. 



rise at large distances. This merely expresses the fact that 
the mass cutoff on the wind size excludes the thinner gas 
that does not contribute significantly to the optical depth of 
the wind. 

Another way of expressing the relation between the ge- 
ometrical cross-section and the efficiency of the winds is by 
the ratio of the emitted La photons to the the number of 
irradiating ionizing photons. This is displayed in figure 8B 
which shows that at small distances, where the winds are 
very dense, the high densities and large optical depths sup- 
press the La. Again, all the models but those of slow and 
dense winds display a similar behavior and the slow and 
dense models are more efficient in converting ionizing pho- 
tons into La and exhibit a rise in this efficiency at large 
distances. 

5.2 The bloated star population 

5.2.1 The integrated line emission 

The line emission was integrated with the cutoff radii 0.1 pc, 
1/3 pc, 1 pc and 56 pc. Figure 9 shows the overall reduced 
X 2 goodness-of-fit score for the integrated line emission with 
an extinction coefficient Eb—v that is required for correcting 
the La /Lip ratio for r ou t = 1/3 pc and 56 pc. The results for 
the O.lpc and lpc cutoffs display the same trends as those of 
the l/3pc with the fit becoming worse as r ou t increases due 
mainly to the over-production of the forbidden lines. The 
extinction coefficient is in the range Eb—v ~ 0.20 to 0.25 
for all the models but the a = 0.5, v* = 5x1 3 cm/s model 
which has the low Eb—v = 0.03. In almost all the wind 
models the reddening correction worsens the \ 2 score, an 
effect due mostly to the increase in the reddened forbidden 
lines. 

Figure 10 A compares the integrated line emission from 
the stellar core to the observed BLR values for the refer- 
ence wind model as function of the cutoff radius. This wind 
model illustrates trends which are common to all the wind 
models: The low Ha, H/? and He I (before reddening cor- 
rection), the decrease of the high density and temperature 
C III A977 line with increasing r ou t and the increase in the 
forbidden line emission with increasing r ou t (The forbidden 
line emission at high cutoffs is outside the graph's scale and 
is comparable to NLR values). The forbidden lines aside, the 
variation of the emission line spectrum with r ou t is not very 
large. Figure 10B shows the integrated line ratios for the 
slow and dense wind model a = 0.5, v* = 5x1 3 cm/s. The 
line ratios, apart of the persistent deficiency of MgIlA2798 
and NvA1240, fit the observed ones quite well. 



increases a larger fraction of the outer part of the wind be- 
comes optically thin. The ratio does not descend smoothly 
with rbh but displays various bumps and local maxima, espe- 
cially at small rbh- This reflects the existence of ionization 
from excited states which can be significant at small dis- 
tances where the gas is highly excited. We note that all but 
the dense and slow wind models exhibit a similar efficiency 
which decreases from ~ 1 at 10 17 cm to ~ 0.1 at 10 20 cm. The 
slow and dense wind models have a consistently higher ac- 
tual to geometric covering factor with some showing a sharp 



5.2.2 Line emission from the innermost BSs 

Figure 11 shows the differential line emission of several lines 
and of the total reprocessed luminosity emitted from a stel- 
lar core containing BSs with the reference model's wind 
structure. The emission exhibits a decrease at small radii 
up to ~ 6xl0 16 cm, then an increase up to 2xl0 18 cm fol- 
lowed by a sharp decline at larger radii. The relatively large 
contribution of the very few innermost BSs can be explained 
by an interplay between the two density regimes of the BS — 
the wind and the star's atmosphere — and the core's stellar 
density distribution function. For optically thick winds the 
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Figure 9. The reduced x 2 goodness-of-fit score and E-g_y. The dark gray column on the right represents the unreddened value, the 
light gray column on the left, the reddened value. The numeric values of the x 2 are given above the columns. E-g_y is given on the left 
of each pair of columns. A) r ou t = 1/3 pc. x 2 based on the BLR values. B) r ou t = 56 pc. x 2 based on the BLR and NLR values. 
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Figure 10. The integrated line emission (relative to La) from the stellar core compared to the observed BLR values. The mean observed 
value is given by the small dot and the range by the error bars. The black triangle, square, pentagon and circle above the error bars give 
the unreddened ratio for cutoffs at 0.1 pc, 1/3 pc, 1 pc and 56 pc respectively. The white shapes beneath the error bar give the ratios 
after being corrected for reddening. A) The reference wind model a = 0.0, = 10 6 cm/s. £/ B _v (0.1 pc) = 0.22, £/ B _v (l/3pc) = 0.20, 
-E B _ v (lpc) = 0.25, _E B _v(56pc) = 0.26. B) The slow and dense wind model a = 0.5, = 5xl0 3 cm/s. £ B _ v (0.1pc) = 0.13, 
-E B _ v (l/3pc) = 0.03, £ B _v(lpc) = 0.08, £ B _ v (56pc) = 0.09. 
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Figure 11. Differential line emission as function of the distance 
from the black hole for BSs with constant velocity wind of v = 
10 6 cm/s. 



total heat reprocessed by a spherical shell of BSs is given by 
the product of the irradiating flux, the wind's cross-section 
and the number of BSs in the shell 

dH(r bh ) ex ^ 2 X7r_R^, x47rrb h ra^(r b h)^r • (26) 



At small radii where R w = Rc < R-k, Rw is almost constant 
(see fig. 6) due to the steep exponential density gradient of 
the star's atmosphere and therefore dH/dr oc ra^(rbh)- For 
the reference wind model this occurs at rbh < 6xl0 16 cm. 
The stellar distribution increases with decreasing radius 
until it reaches a peak at 6xl0 15 cm ~ 2.5 light days. Al- 
though our calculations extend only to as close as 2xl0 16 cm 
from the black hole, we expect that the differential total heat 
emission, following n+, will peak at 2.5 light days and then 
decline sharply. If we assume that the line luminosity ap- 
proximately follows that of the total heat and extrapolate 
linearly to the peak we find that the differential luminosity 
in the lines C ill A977, C IV A1549 and O VI A1035 at 2.5 light 
days is greater than that of all radii up to ~ 10 17 cm, and 
for NVA1240, even up to ~ 10 18 cm. 

The distance where R w = Rc depends on the wind den- 
sity through the parameter M w /v^: The denser the wind, the 
smaller this distance is. At larger radii, where R w > R*, the 
size of a comptonized constant velocity wind such as the 
reference model, is proportional to Nc 1 ^ 2 oc rbh and there- 
fore dH/dr oc 4wr bh n^. The sharp turn-over at 2xl0 18 cm 
reflects the point where begins to fall off faster than r bh . 
The decreasing optical thickness of the wind also contributes 
to the decline in dH/dr. While realistically not all the wind's 
cross-section is optically thick and R w is not strictly con- 
stant at small radii, this qualitative explanation does de- 
scribe the overall behavior of the detailed photoionization 
results. 

We note that unlike the resonance lines and the high 
excitation C III A977 line, the La emission does not increase 
at smaller radii. This is a consequence of the increasing op- 
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tical depth and the line thermalization. It results in large 
line to La ratios in the innermost BSs. 



5.2.3 The number and mass loss of BSs 

Figure 12 shows the fraction of BSs in the core's stellar pop- 
ulation and their number for the various models and the 
l/3pc and 56pc cutoffs. The results reflect the line emission 
efficiency of the individual wind structures. The highest /bs 
values are needed for the fast, highly accelerating wind mod- 
els and the lowest for the slow, decelerating models. These 
trends are also exhibited by the O.lpc and lpc cutoffs. We 
note that the a = 0.0, v* = 5xl0 6 cm/s model, very similar 
to the reference model, requires a high (> 4) fraction of BSs 
at all cutoff radii. 

Figure 13 shows the total mass loss rate and the buildup 
rate of the electron scattering optical depth. The fast and 
accelerated flow models with a < — 1 or v* > 5 x 10 5 cm/s all 
show very high (> lOMg/yr) mass loss rates at all cutoffs. 
At the smallest cutoff, O.lpc, even the densest wind model 
is associated with rather high mass loss and optical depth 
buildup rates of M = 6.5M /yr and f = 0.013 yr" 1 . 

The striking feature in these results is the large differ- 
ence between the properties of the different wind models: At 
a given cutoff radius, /bs extends over a range of 4 orders 
of magnitude, M over a range of 6 orders of magnitude and 
f es over a range of 5 orders of magnitude. 



5.2.4 The BLR size 

Figure 14 shows the flux weighted radius, rf w (H/?), for the 
various models with a cutoff of r ou t = 56 pc. The flux 
weighted radii are much smaller than 56 pc since the stellar 
density falls rapidly beyond 1 pc. The general trends in 
the behavior of rf w (H/?) as a function of the wind model pa- 
rameters a and v+ are also exhibited by all the other lines, 
at all cutoffs. With the exception of the two slow deceler- 
ating models a = 0.5, v* = 5xl0 3 and 5xl0 4 cm/s, there 
is a slow increase in rf w with v+ and a somewhat steeper 
rise with a. This reflects the a and v+ dependence of the 
wind size i? w (n>h) and the line emission efficiency (figures 6 
and 7). The small flux weighted radius of the two slow wind 
models is a consequence of their sharply peaked line emission 
efficiency which reaches its maximum at small radii due to 
the i?mass cutoff on the wind size. rf w (H/?) is, up to a factor 
of 2, also the flux weighted radius of the other lines with the 
exception of the forbidden lines (especially the [O III] A5007 
line) which are emitted at much larger distances. We note 
that the a = 0.5, v* = 5xl0 3 cm/s model does not require 
any cutoff to be consistent with tblr of eq. 20. 

At cutoffs < 1 pc, rf w of the lines is roughly propor- 
tional to the cutoff radius itself. At these small cutoffs the 
difference in rf w of the various lines is small. Depending on 
the wind model and the specific line, the flux weighted radii 
range from a 1/3 to 2/3 of r ou t . 



6 DISCUSSION 

6.1 The scope of this work 

Previous investigations of the bloated stars scenario did not 
follow the path of detailed line emission calculations, but 
rather relied on various simplifying assumptions to predict 
the emission line spectrum (Penston 1988; Kazanas 1989; 
Scoville & Norman 1988). By assuming normal stellar evo- 
lution and the properties of supergiants, or by estimating the 
effects of the irradiating flux on the star, they proceeded to 
check different aspects of the model: The effects of the radi- 
ation on the star, its size and mass loss rate (Mathews 1983; 
Scoville & Norman 1988; Voit & Shull 1988; Kazanas 1989), 
the role of the stars in feeding the black hole and the lumi- 
nosity evolution (Norman & Scoville 1988), the correlation 
between the luminosity and covering factor (Kazanas 1989), 
the line profiles (Scoville & Norman 1988; Kazanas 1989) 
and the implication of the collisional rate (Penston 1988; 
Scoville & Norman 1988; Begelman & Sikura 1991). These 
works did not treat the atomic physics and line emission pro- 
cesses in detail and considered only constant velocity winds, 
with properties similar to those of normal supergiants. In 
particular they did not take account of the varying line emit- 
ting efficiency of the winds and did not employ a detailed, 
self-consistent model of the stellar core. 

Here we adopt a different approach. The AGN envi- 
ronment is drastically different from normal stellar environ- 
ments and its effects on stars and the course of stellar evolu- 
tion are currently unknown. This is to be contrasted with the 
atomic physics involved in the line emission which, although 
computationally complex, are well understood theoretically 
and well tested against observations in various astrophysical 
domains. We therefore use the line emission as the basis of 
our investigation. The use of a photoionization code allows 
us to predict correctly the effects of the ionizing continuum 
on the wind: The comptonization of the outer layers, the 
line emitting efficiency and the emitted line spectrum. The 
known properties of normal supergiants are used only as a 
reference point for a wide range of wind density structures. 
These wind models are then evaluated according to their 
success in reproducing the observed line properties, with no 
a-priori prejudice as to the 'reasonable' wind model. We use 
a stellar core model based on detailed numeric calculations 
(Murphy, Cohn & Durisen 1991). This allows us to assign 
the proper weight to the line emission from different parts of 
the core, to check for consistency with the total number of 
stars in the core and to estimate the magnitude of the col- 
lisional processes once we obtain the number of the bloated 
stars from the photoionization calculations. 

The various components and assumptions of the stellar 
core / bloated star model were discussed in detail in § 2. 
The main input parameters of our models are the choice of 
the stellar core's mass and evolutionary stage, the contin- 
uum luminosity spectrum and the BS's evolutionary mass 
loss rate, wind mass and wind density structure. All other 
physical quantities are derived from these parameters. We 
studied the line emission properties of the single BS as a 
function of the wind structure and distance from the black 
hole with particular emphasis comptonization and line emis- 
sion efficiency. The integrated line emission from all the BSs 
in the core was compared to the observed spectrum. We in- 
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Figure 12. The fraction, . and number, A^gg, of the BSs in the stellar population. The top number at the foot of each column is 
/gg and the bottom one is -/Vgg. A) r ou t = l/3pc. B) r ou t = 56pc. 
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Figure 13. The total mass loss rate from stellar evolution, MS— MS and BS— BS collisions (dark gray columns on the right) and the 
buildup rate of electron scattering optical depth (light gray columns on the left). The numeric values of log(Af) and log(r) are given 
above (if positive) or below (if negative) the respective columns. A) r ou t = 1/3 pc. B) r ou t = 56 pc. 
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Figure 14. Flux weighted radius of the H/5 line. r ou t = 56 pc. 

vestigated the constraints placed on the various wind models 
by the number of BSs and their collisional mass loss rates. In 
addition we studied the possibility of a significant contribu- 
tion from the relatively few BS very near the black hole and 
compared the BLR size, as measured by the flux weighted 
radii of the lines, to that derived from line reverberation 
studies. 

6.2 Main results 

The response of the emission lines to changes of the wind 
structure was shown to be mostly explicable in terms of the 
overall properties of the wind, that is the boundary density 
and the relative sizes of the wind and the star. This inspires 
confidence in our assumption that the results, which are 
based on wind structure models that lack hydrodynamical 
self consistency, are a fair approximation of the line emission 
from realistic bloated stars. The dependence of the emergent 
line spectrum on the conditions at the boundary of the wind 
suggests that the central issue is to understand the mecha- 
nisms that define the outer limit of the line emitting region 
rather than the detailed structure of the wind. 

Broad forbidden lines are not observed in AGN. The 
emission of forbidden lines from BSs, and in particular lines 
of high collisional de-excitation density such as [FeXl] A7892, 
severely restrict possible wind models. Our photoionization 
calculations show that the gas at the comptonized boundary 
has a high ionization parameter (U ~ 3) and emits strong 
broad forbidden lines. Collisional de-excitation becomes ef- 
ficient only very near the black hole, well within the BLR 
radius. This result remains unchanged when a softer con- 



tinuum is used but may be sensitive to the comptonization 
temperature. By artificially lowering the temperature from 
~ 6xl0 4 K to 3xl0 4 the forbidden lines can be suppressed 
up to ~ lpc. Barring that, only mass limited wind models 
whose boundary density is higher than the comptonization 
density avoid excessive broad forbidden line emission. 

The most conspicuous discrepancy between the calcu- 
lated and observed line spectrum is the deficiency of the 
MgIlA2798/La and NvA1240/La ratios, regardless of the 
wind model or BLR size cutoff. (The MgIlA2798/La emis- 
sion of dense wind models does rise up to the marginally 
acceptable value of ~ 0.1, but this emission originates too 
far from the black hole). The MgIlA2798 deficiency can be 
corrected by a combination of a lower optical depth and / or 
a lower ionization parameter. Our photoionization calcula- 
tions assume that the transverse optical depth is infinite 
(slab geometry), which actually over-estimates the true op- 
tical depth of a spherical geometry. 

The NvA1240 deficiency in our models is similar to a 
general trend exhibited by cloud models. Hamann & Fer- 
land (1993) propose that the source of this problem lies in 
the erroneous assumption of a solar chemical composition 
and suggest a composition of higher metalicity for high lu- 
minosity AGN. This is likely to improve the fit but was not 
included in our modeling. The NvA1240 line is strong in 
high density gas (note the steep increase at small radii in 
fig. 5). Its deficiency may be related to our choice of r; n , 
which is perhaps too large and artificially excludes the in- 
nermost line emitting region of the BLR. 

We find, as expected, that the efficiency of the BSs in 
converting the irradiating continuum into line emission de- 
creases with rbh as an increasingly larger fraction of the wind 
becomes optically thin and that this effect becomes more 
pronounced when the wind's density gradient is steeper. 
The efficiency (defined as the ratio between the ionizing and 
the reprocessed luminosities) becomes a sensitive function of 
the mean wind density at high wind densities and reaches 
values of Cf i; 1 for the densest wind model considered 
here. We find that the wind structure, parameterized by the 
wind base velocity and velocity gradient, strongly affects 
EW^(La) by controlling both the wind density and its size. 
The slow decelerating winds have favorable conditions for 
high La emission because the relatively flat density gradi- 
ent results in very large comptonization radii and the high 
density increases the efficiency of the continuum reprocess- 
ing. 

The number of BS required for the BLR is directly re- 
lated to the La emission from the single BS, EW^(La). It 
is therefore also a very sensitive function of the wind den- 
sity structure and changes by more than 4 orders of mag- 
nitude across the parameter space. We find that if we as- 
sume r ou t = 1/3 pc, even the relaxed constraints on the 
fraction of BSs, /bs 1, on the total mass loss from the 
BSs, M is lMg/yr and on the growth rate of the Comp- 
ton optical depth, f es is 0.01 yr -1 limit the acceptable wind 
models to the dense slow models (a > 0, v+ < 0.5km/s or 
a > 0, v* < 1 km/s). Begelman & Sikura (1991) estimate the 
supergiant (SG) collisional mass loss rates for tblr = O.lpc 
and l/3pc at ~ 8500 and ~ SOMg/yr respectively and con- 
clude that SGs cannot play the role of BSs. Our results show 
that these values underestimate in fact the total mass loss 
rate from SG-like BSs (a = 0, v* = 10 km/s) by up to 2 or- 
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ders of magnitude. We therefore also conclude that SGs are 
ruled out by excessive mass loss which is mostly collisional. 

Generally, the number of BS in the successful wind mod- 
els is much smaller than the number of clouds in cloud mod- 
els. The densest wind model (a = 0.5, v+ = 50 m/s) requires 
is 5xl0 4 BSs in the inner 1/3 pc. It is interesting to note 
that Atwood, Baldwin & Carswell (1982) obtained from the 
structure of the Ha and H/? line profiles of Mkn 509 a lower 
limit of 5xl0 4 on the number of individual line emitting 
objects. 

The stellar core considered in this work extends well 
beyond the BLR size that is inferred from line reverberation 
studies. Wind models that produce significant line emission 
at tblr < ^bh < ^nlr require an additional mechanism 
to suppress or obscure this emission. While the need for an 
additional mechanism beyond the scope of the BS scenario is 
not necessarily an objection to this model, we note that the 
dense winds do not encounter this problem since the mass 
cutoff on the size of the wind introduces a natural cutoff to 
the size of the BLR by decreasing Cp(rbh). We find that for 
such wind models this natural cutoff is consistent with tblr- 

Different criteria point to the slow dense models as the 
best candidates for realizing the bloated stars scenario: They 
produce the best fit to the observed emission line spectrum, 
they require a small number of BSs that can be supported 
by the stellar population and avoid rapid excessive mass loss 
and they need no additional mechanism for limiting the size 
of the BLR. The mean wind density (10 8 <; N ^ 10 12 cm/s) 
and its small density gradient resemble those of pressure 
confined clouds and therefore the good fit to the observed 
spectrum is, with hindsight, not surprising. Neither is the 
fact the densest winds are found to be the most efficient in 
converting the ionizing continuum into lines. The similarity 
between the successful wind models and the clouds demon- 
strates the extent to which the observed data and the atomic 
physics constrain any possible BLR models. The most sig- 
nificant result in this feasibility study is that only a small 
number of BSs with dense winds are required for the BLR 
line emission. 

The dense wind models do not resemble solar neigh- 
borhood SGs and to date no such objects were observed 
in any astronomical environment. However, at our current 
level of knowledge, the assumption that the AGN stellar 
population resembles the solar neighborhood population is 
no more likely than the opposite assumption that hitherto 
unobserved bloated stars exist very near the black hole. 



6.3 limitations, uncertainties and future work 

Of the many assumptions that are required to construct a 
physical realization of the bloated stars scenario, we list be- 
low those that we regard as potentially problematic or those 
that were not studied in detail in this work. 

This work is limited to a specific choice of the evolution- 
ary epoch of the core, and with it the black hole's mass and 
the AGN luminosity (§ 2.1). In addition, the MCD stellar 
core model we employ is based on a stellar population that 
does not include BSs. By applying their results to the BSs we 
are implicitly assuming that the dynamical evolution of the 
stellar core proceeds independently of their existence. While 
this may be true for those wind models that require only a 



small fraction of BSs in the stellar population, the existence 
of a significant fraction of BSs must affect the dynamics and 
evolution of the core. 

We have ignored the effects of the irradiating contin- 
uum on the wind density structure and the assumed that the 
wind properties are independent of rbh- The uncertainties in 
the details of the processes that determine the wind bound- 
ary: the choice of the comptonization temperature Tc, the 
amount of mass in the wind, the tidal disruption radius and 
the effect of collisions on the wind, are problematic. The idea 
that mass-limited winds have a sharp, well defined bound- 
ary with a density above the comptonization density plays 
an important role in the collisional suppression of the unob- 
served broad high excitation forbidden lines. It is however 
not clear that such a sharp boundary is physically realistic 
and a more gradual transition from the wind density to the 
HIM density may result in a higher forbidden line emission. 
The actual photoionization calculations also involve several 
assumptions, of which the most problematic is the infinite 
slab approximation. Our model is minimalistic in that it as- 
sumes that the BLR is composed of a single component — the 
bloated stars. We do not consider for example the possibil- 
ity of line emission from dusty clouds or from the accretion 
disk. In addition this model does not address the issue of dif- 
ferent dynamics for different lines (e.g. the observed redshift 
differences between various lines). 

There are two important tests that should be applied 
to the BS scenario in its present level of detail: modeling 
the line profiles and simulating line reverberation mapping 
of the BLR with BSs. The question whether the BS model 
can reproduce the observed profiles and reverberation radii 
is crucial for establishing its viability as an alternative to 
the BLR cloud models. This work is currently under way. In 
addition, steps towards obtaining a hydrodynamically self- 
consistent wind model are required, but achieving physical 
self-consistency in this issue is expected to involve consider- 
able difficulties. 



7 SUMMARY 

We employed a detailed photoionization code and utilized 
a detailed model of the stellar core to study the possibility 
that the BLR clouds are in fact the envelopes of bloated 
stars. The main issues we investigated in this work were 

(i) How is the emission line spectrum affected by the 
structure of the BS? 

(ii) Can bloated stars reproduce the BLR emission line 
spectrum? 

(iii) How many bloated stars are required and can the stel- 
lar population supply them? 

(iv) Can the BSs avoid rapid collisional destruction and 
excessive mass injection to the interstellar medium? 

Our results show that 

(i) The line emission spectrum is mainly determined by 
the conditions at the outer gas layers. This makes 
the physical mechanisms that determine the bound- 
ary density of the wind a central issue, 
(ii) Dense winds (10 8 <; N ^ 10 12 cm" 3 ) can re- 
produce the observed BLR emission line spectrum 
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with a success comparable to that of cloud mod- 
els. The main problems are the under-production of 
MgIlA2798 and Nv A1240. Strong forbidden lines, es- 
pecially [FeXl]A7892, exclude the possibility of lower 
density winds. 

(iii) The required number of BSs is a very sensitive func- 
tion of the density structure of the wind. For dense 
winds it is small enough to be supported by the stellar 
population. BSs with such winds can also avoid rapid 
collisional destruction and the resultant high optical 
depth of the ISM. 
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